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Abstract 

We study the beamline properties of a long baseline neutrino- oscilla- 
tion experiment from CERN to the Pyhasalmi Mine in Finland. We obtain 
the real density profile for this particular neutrino oscillation beamline by 
applying the geophysical data. The effects of the matter density to neu- 
trino oscillations are considered. Also we compare the realistic density 
profile with that acquired from the PREM model. 



1 Introduction 

A future step towards a better understanding of neutrino properties is to build 
long-baseline neutrino-oscillation experiments. One such possibility could be 
to aim a neutrino beam from a proposed CERN Neutrino Factory^] to the 
Pyhasalmi Mine 1 in Finland 2288 km away. The Neutrino Factory and target 
station could be operational after the year 2010. 

The oscillation probability depends on the density of the medium as well as 
the intrinsic mixing parameters. The purpose of this study is to compile a real- 
istic density profile along the CERN-Pyhasalmi neutrino baseline and estimate 
its effect on the neutrino oscillations. 



2 Compilation of the density distribution along 
the baseline 

2.1 The geographical position of the baseline 

The neutrino baseline starts at 46°N 15' 22.93", 06°E 03' 01.10", +450 m above 
sea level, and ends at 63°N 39' 34.74", 26°E 02' 29.94", -1345 m. Its total 
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Figure 1: Location of the CERN-Pyhasalmi baseline (solid black line) on the 
map of main tectonic elements of Western Europe^]. The location of the EGT 
in the swathe between the two broad lines. 



length is 2288 km. The geographical position of the baseline and its penetration 
depth were calculated using the geocentric Cartesian coordinate system and 
then transformed to the geographical coordinates and depth with respect to the 
WGS-84 ellipsoid [2] using the transformation equations in !3j. The location of 
the baseline is shown on the simplified tectonic map of Europe (Fig. The 
baseline goes mainly through the lithosphere (the Earth's crust and the upper- 
most part of the Earth's mantle) and crosses the main geological structures and 
boundaries of the Europe. Its deepest point (103.81 km) is located below the so- 
called Trans-European Suture Zone (TESZ). Thus, the density variations along 
the line can be due to different thickness of the crust and different density of the 
crust and lithospheric mantle within various tectonic units. Another affecting 
factor is the depth to the boundary that separates the non-convecting litho- 
spheric mantle from partially molten convecting and less dense asthenosphere. 
This information can be obtained from results of recent lithospheric studies in 
Europe. 

As seen from Fig. the large part of the baseline is located within the 
study area of the European Geotraverse project (EGT). The EGT was a 4600 
km long and 200 — 300 km wide lithospheric transect across Europe from Nor- 
way to Tunisia. The multinational, multidisciplinary research resulted in a 
comprehensive cross-section of the European lithosphere to the depth of 450 
km . It revealed a significant contrast between the thickness and structure of 
the lithosphere of younger western Europe (less than 90 km) and the old cold 
lithosphere of Fennoscandia and eastern Europe (more than 160 km), occur- 
ring beneath the TESZ (Fig. EJ. Later, the more detailed information about 
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Figure 2: A sketch demonstrating the location of the CERN-Pyhasalmi baseline 
on the composite cross-section of the lithopshere along the EGT jl| (solid black 
line). The depth to the lithosphere-asthenosphere boundary beneath the Trans- 
European Suture Zone is modified in accordance with the S-wave velocity model 
in 0. M, Moho boundary; TZ Tornquist Zone; TEF, Trans-European Fault. 



the lithosphere-asthenosphere boundary beneath the TESZ was obtained within 
the seismological EUROPROBE/TOR project 00. 

In Finland, the CERN-Pyhasalmi baseline is located within the study area 
of the EUROPROBE /S VEK AL APKO project 0. As a part of the project, 
the 3D density model of the SVEKALAPKO area on a 2 km x 2 km x 2 km 
rectangular grid down to the depth of 70 km was obtained ■ 

Thus, the abundant geophysical information about structure of the crust 
and upper mantle obtained within three large geoscientific projects in Europe 
mentioned above (e.g. EGT, TOR, SVEKALAPKO) allows us to compile a 
realistic density profile for the CERN-Pyhasalmi baseline. 

2.2 Composition and density of the continental crust and 
the lithospheric mantle 

The thickness of the crust in Europe varies from 26 km in its western part to 65 
km beneath the Fennoscandia. The most abundant elements in the continental 
crust are O (46.4 wt%), Si (28.15 wt%), Al (8.23 wt%), Fe (5.63 wt%), Mg 
(2.33 wt%), Ca (4.15 wt%), Na (2.36 wt%) and K (2.09 wt%) 0, occuring in 
the Earth in the form of various silicate minerals. The density in the continental 
crust generally increases with depth due to the lithostatic pressure resulting in 
rock compaction and also due to decrease of the content of Si02 in the rock- 
forming minerals of lower crustal rocks. The strongest density contrasts in the 
crust exist between the sedimentary cover and the bedrock and also at the so- 
called Mochorovichich boundary (Moho boundary) separating the crust from 
the mantle. 

The major part of the baseline is located within the upper 100 km of the 
Earth's litospheric mantle. The physical properties and composition of it are 
known from geophysical studies and from direct measurements on samples of 
mantle rocks that have been overthrusted or exhumed to the surface by vari- 
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ous tectonic and magmatic processes. The main elements in the upper mantle 
are Fe, Mg, Si and 0. The main rock forming minerals for upper mantle rocks 
(peridotites) are olivine ((Mg,Fe)2Si02), orthopyroxene ((Mg,Fe)Si02), clinopy- 
roxene ((Ca,Na)(Mg,Fe,Al, Cr,Ti)(Si,Al) 2 6 ) and spinel (MgAl 2 4 ). It is be- 
lieved that the continental lithospheric mantle has undergone significant melting 
through geological time and is depleted in such components as Fe, Al, Ca and Ti. 
As a result, the density of the continental lithospheric mantle generally depends 
on the tectonothermal age: younger lithospheric mantle is less depleted and 
hence is denser. The conversion from fertile to depleted mantle is expressed by 
a decrease in clinopyroxene and orthopyroxene and relative increase in olivine, 
and also by an increase of Mg content. That is why the lithospheric mantle 
beneath the young Phanerosoic Western Europe has the average density of 3.37 
g/cm 3 , while the older (Proterozoic) mantle beneath the Fennoscandian Schield 
is less dense (3.34 g/cm 3 ) [To] . 

In comparison to the lithospheric mantle, the asthenosphere is fertile and 
chemically more homogeneous because of convection. The density of the as- 
thenosphere is less than that of the lithospheric mantle and depends mainly on 
the density and content of partially molten material (tholeitic basalt) Ql] . 

The amount of partial melt in the asthenosphere can be estimated from the 
stability condition for the olivine-molten basalt mixtures. They are mechanically 
stable only in the case when the melt is concentrated within isolated, non- 
connected inclusions (pockets) . The theoretical modelling of elastic and electrical 
properties of such mixtures^2] demonstrated that in some cases only 5% of melt 
inclusions is enough to form a perfectly interconnected network. This value (less 
than 5%) is in agreement with estimates obtained by teleseismic tomography 
studies ■ Thus, if the melt content is less than 5%, and the density of 
the molten basalt under upper mantle pressure-temperature conditions is 2.72 
g/cm 3 [2], then the lithosphere-asthenosphere density contrast is less than 0.04 
g/cm 3 . 

The lithosphere-asthenosphere density contrast beneath the TESZ can be 
roughly estimated also from the S-wave velocity (V s ) model in 0. In accordance 
with it, the V s in the mantle lithosphere is 4.57 km/s, while the V s in the 
asthenosphere is 4.36 km/s. Using the scaling factor relating decrease in V s to 
the decrease in density, that is equal to 0.05 at a depth of about 100 km |15j . 
the asthenophere density beneath the TESZ is 3.34 g/cm 3 . 

These estimates of the asthenosphere density agree with the results of re- 
gional (medium-wavelength and long-wavelength) gravity studies in Europe 
(3.2 — 3.33 g/cm 3 ) JUl ■ Summarising all the data about the density in the 
asthenosphere cited above, we can conclude that the range of possible values of 
the asthenosphere density for Western Europe is 3.2 — 3.34 g/cm 3 . 

2.3 Geological setting and density values along the base- 
line 

From to 175 km the baseline goes through the Earth crust beneath the Jura 
Mountains, Molasse Basin and the Upper Rhine Graben. The thickness of the 
crust here is 30 — 28 km JZj- The first 21 km of the line are located within 
the sedimentary cover that is 3300 — 3400 m thick^H]. The density within 
the sedimentary cover composed of limestones, sandstones and shales gradually 
increases with depth from 2.4 to 2.6 g/cm 3 . 
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The next portion of the line (22 — 175 km) goes through the upper(22 — 71 
km), middle (72 — 119 km) and lower crystalline crust (120 — 175 km) with the 
densities of 2.7-2.73 g/cm 3 , 2.8-2.87 g/cm 3 and 2.96-2.97 g/cm 3 , respectively. 
The thickness of the layers and their density were taken from the 3-D gravity 
studies [El (201 ■ 

The next part of the baseline (175 — 732 km) goes through the Palaeosoic 
lithospheric mantle (3.37 g/cm 3 ). Basing on studies by I20jj21j we assumed 
no upwelling asthenosphere beneath the Upper Rhine Graben (175 — 430 km). 
Thus, the baseline goes through the uprising asthenosphere only at 732 — 1019 
km. As it was shown in previous Chapter, the density of the asthenosphere 
here can be 3.20 — 3.34 g/cm 3 . The lithosphere-asthenosphere boundary along 
the baseline was estimated in accordance with the lithosphere thickness map[22 
and from the model|3, with the accuracy of 50 km. 

Prom 1020 to 1956 km the baseline is located within the Proterozoic litho- 
spheric mantle (3.34 g/cm 3 ), then it returns to the crust at a depth of 52.01 km. 
The density values along the final part of the baseline (1957 — 2288 km) were 
taken from the 3-D density model of the SVEKALAPKO area[H]. The crust in 
this part of Finland consists of four major layers: the upper crust(2.73 — 2.81 
g/cm 3 ), the middle crust(2.86 - 2.89 g/cm 3 ), the lower crust(2.86 - 2.89 g/cm 3 ) 
and so-called high-velocity lower crust(3.01 — 3.2 g/cm 3 ). The density values 
along the baseline are summarised in Figure 21 
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3 Neutrino matter oscillations along the baseline 

Next we will discuss the effect of the realistic CERN-Pyhasalmi density pro- 
file to neutrino-oscillations. In a general three flavour framework the distance 
evolution equation of the neutrino flavour states can be written as [23] 

i Tx j = \h UM2u] + 1 1 ) >< i ^ i (i) 

and we use the PDG endorsed parametrization_24| for the unitary mixing ma- 
trix U. The potential term AV = ±^/2GpY e p arises from the coherent for- 
ward charged-current scatterings of electron neutrinos from the matter elec- 
trons. Throughout this paper we assume Y e — 0.5, because of the crust and 
mantle material content, which is mainly Si and O. We solve the equation Q 
numerically by using a fourth order Runge-Kutta -algorithm. 

At energies relevant to a neutrino factory beam (above 1 GeV) , the v e — > 
-oscillations are governed mainly by the mixing angle #13 and the mass squared 
difference 8m%^ ~ 5m\ 2 . For typical values of Sm^ = +3.0 ■ 10~ 3 eV 2 , #13 = 
0.1 and density p w 3.3 g cm -3 , the matter resonance energy in the leading 
approximation's] is about ~ 11 GeV, but the distance 2288 km is much smaller 
than the corresponding oscillation lenght (L osc ~ 45000 km) and therefore no 
total resonance conversion can happen. Our simulations show that the muon 
neutrino appearance probability at the first oscillation maximum E ~ 5 GeV 
is enhanced by about a factor of two compared to the appearance probability 
in vacuum. For a negative sign of 5m\i, corresponding to the inverted mass 
hierarchy, the appearance probability is suppressed by about a factor of two. 
For V e — > F^-oscillations the situation is opposite (shown in Figure EJ. 
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Figure 3: An example of muon neutrino and antineutrino appearance probabil- 
ities for different energies and mass hierarchies. (Parameters set to Sm^i = 
5 • l(T 5 eV 2 , <5m| 2 = ±3.15 • lCT 3 eV 2 , sin 2 26»i2 = 0.87, sin 2 26> 23 = 1.0, 
sin 2 20i3 = 0.1, <5 = 0) 



The largest uncertainty in the obtained density profile is the astenosphere 
at L = 733 . . . 1010 km, whose size and density can vary by AL — ±100 km and 
p = 3.2 . . . 3.34 g cm~ 3 . Because of this uncertainty, we consider three possible 
density profiles (shown in Figure 01 together with the profile acquired from the 
PREM-model[2Sl): an average density and average astenosphere width profile 
with p as t = 3.27 g cm -3 (1), a wide astenosphere with low density p as t = 3.20 
g cm~ 3 (2) and a narrow astenosphere with high density p as t = 3.34 g cm -3 
(3). We compare the muon neutrino appearance probabilities calculated with 
the different matter profiles. The absolute differences are shown in Figure 
We see that, for an upper limit value[22j sin 2 2(9i3 = 0.1, the differences between 
the appearance probabilites at oscillation maximae are of the order of few 10~ 4 . 
This corresponds to a relative error below 1%. We conclude that the errors of 
the uncertain nature of the astenosphere can be considered as a small, often 
negligible, error in future analysis. 

4 Discussion 

The PREM density profile is commonly used in Earth matter density neutrino 
oscillation studies. There has been discussion j2Hj about the accuracy of this ap- 
proximate model. Therefore we also study the oscillation probability differences 
between the PREM and the realistic density profile. 

We consider the v e — > -oscillations. The absolute difference for the muon 
neutrino appearance probability between the two density profiles is shown in 
Figure El The difference of the appearance probability between the two density 
profiles is, around the first oscillation maximum, about a factor of 3 times 
larger than the error due to the astenospheric uncertainties. So, in obtaining 
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Figure 4: The different density profiles used for error analysis. 
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Figure 5: The effect of the astenospheric uncertainties to the muon neutrino 
appearance probability. Absolute differences of matter profiles 2 and 3 from the 
density profile 1 as a function of energy. (Parameters set as in Figure with 

Sm 2 32 > 0.) 
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Figure 6: The absolute difference of the muon neutrino appearance probabilities 
P(PREM) - P(REAL). (Parameters set as in Figure 0) 



the realistic density profile, the errors of neutrino oscillation simulations are a 
small step better, allthough they are not drastically improved. 

Random matter fluctuations of the density profile have been considered in 
Ref. (221 ■ It is estimated that with different realistical random matter fluctu- 
ations the average error of P Me for a distance L w 2300, average density p = 3 
g cm" 3 and energy E v = 30 GeV is |AP Me | ~ 0-5 . . . 1 • 10 . In our case (with 
parameters as in Ref. .2Q|), the difference between the PREM and the realistic 
case is |AP Me | = 0.4 • 10~ 5 and lAP^I = 0.6 • 10~ 5 . Hence it seems that the 
error result from the random matter fluctuation model fits nicely in to our 2288 
km baseline. 

We will also like to point out that the calculated relative error between the 
neutrino factory beam muon event rates (from v e — > v^) using the PREM and 
realistic density profiles is below ~ 2% for beam energies around 1 — 10 GeV 
and is considerably smaller for larger energies. 

We conclude that we have made up a reference model for the CERN- Py- 
hasalmi baseline that is accurate enough for all studies of long baseline neutrino 
oscillations. Although the results are quite similar to the PREM model, there is 
no excuse for using a less accurate model. However, for other baselines that have 
not been specifially modelled the PREM model may be a good approximation 
to start with. 

The calculated event rates and parameter reaches of the CERN-Pyhasalmi 
neutrino beam will be discussed in a future paper. 
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